Dispersion interaction between two atoms in electromagnetic fields 
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We present a new theory of atom-atom dispersion interaction in the presence of electromagnetic 
fields. The theory takes into account the absorption and emission of virtual photons leading to the 
resonance contributions to the interaction potential in the case of non-equilibrium dynamics. 
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Recent progress in experimental methods a 1 1 a 
enabling one to measure the thermal component of the 
Casimir-Lifshitz force [3], stimulated a new wave of the- 
oretical works in Casimir, Casimir-Polder (CP), and van 
der Waals dispersion forces at the presence of electro- 
magnetic field [3, 0, 13, 13] ■ The problem of the Casimir- 
Lifshitz interaction between real metals is still not re- 
solved. Although, the interaction is described by the 
famous Lifshitz theory, the theoretical predictions dra- 
matically depend of the model describing the real metal 
(see for example 13, II , 12 1 and references therein) . An- 
other problem to be clarified is the non-equilibrium CP 
interaction between a single atom and a metal (or dielec- 
tric) plate if the temperature of the plate differs from the 
temperature of the atomic gas. This force was recently 
measured with the help of spectroscopy of an atom in- 
teracting with the dielectric slub [13], but the agreement 
with the theoretical predictions has not been achieved 

0- 

The theoretical treatment of the CP interaction under 
non-equilibrium conditions is based either on the linear- 
responce (or Lifshitz) theory 0, H, Q3] , or macroscopic 
quantum electrodynamics [l5|. As it has been demon- 
strated by Buhmann and Scheel [13] one should distin- 
guish two different problems. If the atom is thermalized 
and coupled to its thermal bath, the both approaches 
lead to the same result. If the atom is not thermalized, 
the results of the linear-responce theory could underes- 
timate the the CP force [l6|. The macroscopic QED 
approach takes into account the possible absorption and 
spontaneous/stimulated emission of a thermal photon by 
the interacting atom. These contributions result in the 
resonance enhancement of the CP force even if the atom 
is in its ground state. For case of thermalized atom, this 
terms cancel out and one arrives at the predictions of 
Lifshitz theory fvfy . 

To understand the physical mechanism of the CP inter- 
action under non-equilibrium conditions it is important 
to have a clear picture of the simplest case - the dispersion 
interaction between two atoms in electromagnetic fields. 
Milonni and Smith treated such kind of interaction [13] 
with the help of source theory [ljj. They showed that the 
interaction potential between two atoms in electromag- 
netic field could be obtained from the vacuum field poten- 
tial by the simple substitution Uuik/2 hu^ (A^a + 1/2) 
[13], where Wk is the energy of the photon of {k, A} 
state, with k and A the wave vector and the polariza- 



tion index of EM field, N k \ is the number of photons 
in the state {kA}. Deriving their formula, Milonni and 
Smith took into account the scattering of electromagnetic 
field on atoms but they ignored the virtual absorption 
and/or spontaneous/stimulated emission of photons by 
the atoms. Obviously, this approach is valid for the case 
of equilibrium due to the detailed balance between the ra- 
diation and the state of atoms. But if the atoms and the 
radiation are not in the equilibrium their result is in con- 
tradiction with the recent theoretical results obtained for 
CP force between an atom and a dielectric slub 13, 13] 
taking into account virtual absorption and emission of 
photons. 

Recent theoretical and experimental progress in non- 
equilibrium Casimir physics inspired us to reexamine the 
interaction between two atoms in electromagnetic field. 
To calculate the interaction potential we use a QED 
method based on the Keldysh Green function technique 
[i3 |. We start with the derivation of the general formula 
for the interaction between a ground state atom and a 
dielectric medium in electromagnetic field. To simplify 
the calculations we suppose that the number of photons 
of the external EM field is negligible at the transition 
energies of the ground-state atom. This simplification 
enables us to treat the atom as a ground-state one dur- 
ing the interaction and calculate the interaction potential 
as the energy shift of the atomic level. Then we consider 
the interaction of a ground-state atom with a medium 
at non-zero temperatures, provided the atom is not ther- 
malized. We rederived the Buhmann-Scheel formula [13]. 
Next, we consider the interaction between two atoms in 
electromagnetic field and demonstrate that the absorp- 
tion or spontaneous emission of photons play major role 
for non-equilibrium case resulting in resonance terms of 
the interaction potential. We show that the dispersion 
force for an excited atom and a ground-state one is en- 
hanced by factor N(ujb) + 1, where lub is the transition 
frequency for the excited atom, compared to the Power 
and Thirunamachandran result for electromagnetic vac- 
uum [ill. [22!] . If both atoms are in their ground states, 
the difference between the vacuum case and the case of 
external field is more dramatic. The presence of the elec- 
tromagnetic field results in the resonance force, which 
drops like R~ 2 with the distance between the atoms R 
at the retarded regime (R >> 2ir/k). We discuss the 
validity of the Milonni-Smith formula. 

We start with the CP interaction between a single 
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ground-state atom (atom A) at a position R^ and an 
arbitrary dielectric medium at the presence of electro- 
magnetic field. If the number of photons at the transition 
frequencies of the atom is negligible, the atom does not 
change its internal state. In this case we can treat the 
CP potential as the energy shift of the ground level of 
the atom U = AeA [HI- The Hamiltonian of the system 
atom A - electromagnetic field - medium reads 



H = Yj eib l bi +J2 k ( a kA"kA + 1/2) + H med 



kA 



V> + (X)dE(X)iP(X)dr 



(1) 



Where ej is the bare energy of atom A at state i, hi 
is the annihilation operator of this state, k is the wave 
vector of electromagnetic field, A is the polarization in- 
dex, akA is the annihilation operator of state {k, A} of 
the electromagnetic field, H m ed is the Hamiltonian of 
the medium interacting with the electromagnetic field, 
"4>(X) = J2 i 4>i{r, Ra)&i with <f> wavefunction of the atom 
A, X = {r, t}, E is the operator of electromagnetic field, 
d is the dipole moment of the electron. 

Following the earlier treatment for a dipole interaction 
of a ground state atom with a dielectric medium in elec- 
tromagnetic vacuum we generalize our previous re- 
sult to the case of electromagnetic field surrounding the 
atom and the medium The retarded Green func- 

tion of atom A obeys the Dyson equation. In derivation 
we followed the standard Keltysh method for the electro- 
magnetic field and a modified one for the atoms 2(], 2~3 |. 
Treating the atom, we directly implemented Wick's the- 
orem for a single atom making use of the fact that the 
normal products of the atomic operators of all orders but 
the second one are zero. The normal product of the sec- 
ond order is just the density matrix of the atom A 
In energy representation the Dyson equation reads 

G r (w,r,r / ) = G?(w,r,r / ) 

+ J G° r (u;,r,r 1 )M r {u),r 1 ,r 2 )G r (cj,r 2 ,r')dT 1 dr2 (2) 

where G° is the retarded Green function for a non- 

(3) 



interacting atom, M r is the mass operator [24 1 

M r {X, X') = id v d v 'G r {X, X')D VU ' {X\ X) 
Where D is the casual photon Green function 

D UU '(X,X') = -i{fE u '{X')E»{X)) 



(4) 



T is the time-ordering operator Neglecting the transition 
of the atom A to its exciting state and using the pole 
approximation we arrive at the solution of the Dyson 
equation |(2]) 



(5) 



with Mf = J (j)*(r)M r (v,v')(j) i {v')dvdv' . Thus, the en- 
ergy shift of the ground level of atom A is Ua = ReM". 
Following the calculations of [13, US], we obtain 



U A 



-Re- 



2tt 



<x v a v {u)D vv (w,RA,R A )dw (6) 



This formula generalizes the one [20j] obtained for the in- 
teraction of an atom with a dielectric medium at zero 
temperature if the electromagnetic field is in its vacuum 
state. Now we include the external EM field. The polar- 
izability of the atom at the internal state k is given by 
the standard formula (we put h = l,c = 1 throughout 
the paper) 



V V t \ 



UJ 



iO uijk + u) + iO 



(7) 



with ujkj transition frequency of the atom from state k to 
state j, dj k the ^-th projection of electric dipole matrix 
element between states j and k, iO describes the analyt- 
ical properties of the polarizability. 

To calculate the Green tensor we use the standard 
Keldysh technique 25, 2(|. First we notice that the ca- 



sual Green tensor can be written as a sum of the retarded 



Green tensor and the Green tensor D12 25j, [2 



D = D r + D u 



where 



D^{X,X') = -i{E»\x')E»{X)) 



(8) 



(9) 



Substituting © and © into © we find 



U A = -Re— J a A »(cj)D™ (co,R A ,R A )du 

if 00 , 

- Re 2^J Q oT A v (u)D% (u>,R A ,R A )du (10) 



The first term of the Eq. lfTOj) is the standard interaction 
potential for the ground-state atom interacting with the 
vacuum electromagnetic field at the presence of an arbi- 
trary dielectric body [27j. The second term describes the 
interaction of the atom with the photons of the exter- 
nal electromagnetic field or the photons radiated by the 
media. 

For zero temperature the density matrix D12 describes 
only spontaneous emission by the medium and Ex. iflOj) 
results in Ex. (41) of 0. 

As the first example we consider a ground-state atom 
embedded in the system electromagnetic field - dielec- 
tric medium at thermal equilibrium at temperature T. 
We suppose that the atom does not change its initial 
state, it means that the time scales are short compared 
to the inverse ground-state heating rates of the atom F^ 1 . 
This situation was recently considered by Bushmann and 
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Scheel fl6| . The casual Green tensor of the thermal elec- 
tromagnetic field is [28| 

D{uj, R A , R A ) = ReD r (co, R A , R A ) 
+icoth(uj/2T)ImD r {u},'R,A 1 ~RA) (11) 



It can be rewritten as 



D = (2N(u>) + l)D r - 2N(uj)ReD r 



(12) 



where N(w) is the occupation number of the photons at 
frequency uj Using Ex. J6]), we find for an isotopic atom 

oo 

U A = T^2(1- -5 m0 )aY{iU)D v r v (if m ) 

ro=0 

-l/a^JV^IdypiteDrCWjfcJ^-fc) (13) 



where £ m = 2irmT is the Mazubara frequency. This re- 
sult coincides with the Ex. (25) of [lfj]. The first term 
describe the standard Lifshitz result for the interaction 
of a thermalized atom with the medium at thermal equi- 
librium. The second term, as it was pointed out by Buh- 
mann and Scheel, is the resonance contribution to the 
Casimir-Polder force due to absorption of thermal pho- 
tons by the atom, which is not thermalized 

Let a ground-state atom A be at a position R A and 
atom B at a position Rg. Atom B can be either in ex- 
cited or ground state. For simplicity we suppose that 
the atoms are exposed to isotropic unpolarized electro- 
magnetic field, i.e. depends only on |k| = uj , 
Nkx = N(lj). We suppose that N(uf k ) = to simplify 
the calculation. We will attack the problem perturba- 
tively. The density matrix Di 2 can be calculated with 

According to ([6]), 



25, 2 



the help of Keldysh technique 
we should integrate the Green tensor with respect to pos- 
itive frequencies u> > 

D 12 (X,X')=D° 12 
+ J D° r (X,X 1 )IL r (X 1 ,X 2 )D° 2 (X 2 ,X')dX 1 dX 2 

+ J D° 2 (X,X 1 )Tl a (X 1 ,X 2 )D a (X 2 ,X')dX 1 dX 2 

-J D° r {X,X 1 )IL 12 (X 1 ,X 2 )D° a (X 2 ,X')dX 1 dX 2 (14) 

Here we suppressed the tensor indices v. D° means free 
photon Green tensor. X = {r, t} The polarization oper- 
ators II describe interaction of electro mag netic field with 
atom B which is in its initial state n [20j 

TLf{u,r,r') = -a v /(oj)5(r-r') = C'*( W ,r',r) 
11%' (oj, r, r') = ]T 2nid v n ld» m *5(cj - io nm )8(v - r')(15) 

m 

The polarization operators fl r and n a describe elastic 
scattering of light on the atom B, while II12 describes 



the downward transition of the atom as a result of spon- 
taneous or stimulated emission. For positive frequencies 
(uj > 0),the free field Green tensor D\ 2 we can repre- 
sented as 

D° 12 (u>, r, r') = N(uj)(D°(uj, r, r') - D° a (co, r, r')) (16) 

The formula JH]) along with {T5j) and (JTeJ) yield 

Dia(w, Ra, Ra) = D° 12 {u, R A , Ra) 
+2N(iu)D r (oj, R A , R B )U r (uj)D° r (^, Rb, Ra) 
-2N{u>)ReD Q r (w, R A , R B )II r (a;)D°(a;, Rb, Ra) 
+N(cj)D° r (ij, R A , R fl )n 21 ( W )D°( W , R B , R a ) 
-(N(oj) + l)D° r (uj, R A , R B )U 12 (uj)D a r (uj, R B , R A ) (17) 

The polarization operator II21 describes the absorption 
of a photon by atom B and an upward transition of the 
atom B. 



n 2 T>,r,r') 



£ 27riOO(w ~ "rnn)S(T ~ r')(18) 



The first term of the r.h.s of Eq. (fT7]l describes the free 
electromagnetic field. This term results in the optical 
Stark shift of the energy levels of atom A and we omit 
this term. The second and the third terms describe the 
elastic scattering of light on atom B. The fourth term is 
responsible for absorption of a photon by atom B and the 
last term is for the spontaneous and stimulated radiation 
of a photon by atom B. For isotopic atom B prepared 
in a mixed state with the probability to find the atom in 
n-th state p n we obtain 

/>oo 

U a = Re— I duj(N(u) + 1/2) a A (w)a B (w) 

T JO 

x(^(co,R A R B )) 2 
+ - Re J2 I df nn 1 2 a A (w m „ ) | D° r (w mn , R A , R B ) | 2 

mn 

x [p n N(uj mn )8(uj mn ) - p n (N(oj nm ) + 1) 6(u nm )] (19) 



The first term of the Eq. I|T9|) coincides with the one ob- 
tained by Milonni and and Smith for interaction between 
two atoms in electromagnetic field using the source the- 
ory [l3|. The second term describes the absorption of a 
photon by atom B (first term in the square brackets) as 
well as spontaneous and stimulated emission by atom B 
(second term in the square brackets). The second term 
is resonant. It could be easily checked by substitution 
of the polarizability of atom A (Ex.jl])). Let us consider 
a case of thermal equilibrium. The probability to find 
an atom in its n-th state p n is given by the Boltzmann 
distribution, while the number of photons obey the Bose 
distribution. One can easily check that in equilibrium, 
the second term is zero due to detailed balance, and we 
find that the equation (fl9]) coincides with the Lifshitz 
formula fl7ll. 
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Thus we arrive at a conclusion that the Milonni for- 
mula describes only the interaction between two atoms 
at equilibrium. But for the non-equilibrium case, if one 
or both atoms are not thermalized, the Milonni formula, 
which does not take into account the absorption and 
emission of photons, underestimates the interaction po- 
tential between two atoms in the electromagnetic field 
even if the atoms are in their ground states. 

Now to obtain quantitative results, we substitute the 
explicite expression for the retarded Green tensor [28] 



XLU 



I 

ujR 



M) 5 



3i 
ujR 



- 1 



exp(iujR) 
R : 



(20) 



where R = |r — r'\ and s u = (r — r') u /|r — r'|, into (fl9|l 
For the non-resonance and the resonance parts of the 
potential we find 

1 f°° ( 2 

U An r = 57 / a A {iu)a B {iu)u A 1 + — 

ttR Jq \ uR 

5 6 3 



{uRf (uRf {uR) 4 i 
x (2N(iu) + 1) eyq>[-2iuR]du (21) 

4 \ d kj\ 2 \ d mn\ 2uJ kj(A)^tn(B) 

= 9R 2 ' 



U A 



i ,2 _ , ,2 

kj(A) W m»(B) 



X 1 + 



1 



i^kj(A)R) 2 {Ukj(A)R) A 
X [pnN((J mn (B))0(Umn(B)) 
X -p„(N(u) nm ( B )) + l)0{Unm(B))] 



(22) 



If the electromagnetic field is in its vacuum state 
(N(lo) = 0)) the first term in the square brackets of (|24|) 
is zero. The second term is not equal to zero only if atom 
B is excited. In this case we obtain the result by Power 
and Thirunamachandran 12 ll. l22l . 



U Ares 



4 ^K/l 
9R?2-^ 



J kj(A) 



1 



{ukj(A)R) 2 {vkj(A)RY 



mn(B) 
#(w n m(B))-(23) 



For the limit R >> 2n/k, when the retardation effects 
are significant, the potential drops with the distance be- 
tween the atoms as Uatcs oc R~ 2 (compare to the non- 
resonance contribution which is given by Casimir-Polder 
formula U Anr oc R~ 7 ), this weird behavior of the po- 
tential was discussed by Power and Thirunamachandran 



21|,|22|. (see also [20j, where the interaction potential be- 
tween an excited atom and a ground-state one embedded 
in a dielectric medium was calculated for vacuum state 



of electromagnetic field) . In the presence of electromag- 
netic field, the resonance contribution to the interaction 
between a ground-state atom and an excited one is en- 
hanced by factor N(cu B ) + 1 compared to the interaction 
potential in electromagnetic vacuum. 

If both atoms are in their ground state, we should take 
into account only the first term in square brackets of f2~4|) 
For the retarded regime (R » 2n/k) 



UAres 



4 ^K?\d 
9R 2 ^ 



nij kj(A) 
XN(uJ mn (B))d(w mn (B))- 



7U7l(B) 



(24) 



This potential is resonant and drops as R~ 2 with the dis- 
tance. It means that the contribution of the resonance 
term due to photon absorption to the interaction poten- 
tial in non-equilibrium situation could be much greater 
then the contribution of a non-resonance one described 
by the Milonni formula. We should mention, that the 
approach we developed in this presentation is valid for 
the initial stage of interaction when the atoms and elec- 
tromagnetic field are not in the equilibrium. Obviously, 
after the equilibration the resonance contribution to the 
interaction potential is zero and we obtain the Milonni 
formula. 

In conclusion, we reexamined the atom-atom disper- 
sion interaction in electromagnetic fields. We showed 
that the absorption and emission of virtual photons 
results in the resonance contributions to the non- 
equilibrium interaction potential. 
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